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Abstract 



We determine the weak coupling \V c i\ between the b and c quarks using a 
sample of 3 million BB events in the CLEO detector at the Cornell Electron 
Storage Ring. We determine the yield of reconstructed B — > D*£v decays as 
a function of w, the boost of the D* in the B rest frame, and from this we 
obtain the differential decay rate dF/dw. By extrapolating dT/dw to w = 1, 
the kinematic end-point at which the D* is at rest relative to the B, we 
extract the product |Vd,|.F(l), where ^"(1) is the form factor at uu = 1. We 
find |L c fe|^(l) = 0.0431 ± 0.0013(stat.) ± 0.0018(syst.). 
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Typeset using REVTeX 



In the Standard Model, the weak decay of quarks is described by a unitary 3x3 matrix 
HJ. This CKM matrix describes the flavor mixing among the quarks and the amount of CP 
violation through its single non-trivial phase. Precise determinations of the CKM matrix 
elements are essential to test the Standard Model. This Letter presents an improved mea- 
surement of \Vcb\, the coupling of the b quark to the c quark. The CKM matrix element \V c b\ 
sets the length of the base of the unitarity triangle (UT), which displays one CKM unitarity 
condition, and normalizes the constraint on the UT from indirect CP violation in K° decay. 

One strategy for determining \V c b\ uses the decays B° — > D* + £~v and B~ — > D*°£~v. 
The rate for these decays, however, depends not only on \V c b\ and well-known weak decay 
physics, but also on strong interaction effects, parameterized by form factors. In general, 
these effects are difficult to quantify, but Heavy Quark Effective Theory (HQET) offers a 
method for calculating them at the kinematic point at which the final state D* is at rest with 
respect to the initial B meson. In this analysis || , we take advantage of this information: 
we divide the reconstructed candidates into bins of w, where w is the scalar product of the B 
and D* four- velocities. Using these yields, we measure the differential rate dT/dw for w > 1, 
and extrapolate to obtain the rate at w = 1, which, combined with theoretical results, gives 
\V cb \. 

We analyze 3.33 million BB events (3.1 fb ) produced on the T(45) resonance at 
the Cornell Electron Storage Ring (CESR) and detected in the CLEO II detector ||. In 
addition, we use 1.6 fb _1 of data collected slightly below the T(4S') resonance for the purpose 
of determining continuum e + e~ — > qq(q = u, d, s, c) backgrounds. 

To identify D* candidates, we reconstruct the decay chains D* + — > D tt + and D*° — > 
D°7T followed by D° — > K~tt + . We first combine kaon and pion candidates in hadronic 
events to form D° candidates. Signal candidates lie in the mass window \m{K^) — 1865 1 < 
20 MeV/c 2 . We then add a slow vr to the D° candidate to form a D*. For D* + (D*°) 
candidates we require Am = m(K7nr) — m(Kit) to be within 2 MeV/c 2 (3 MeV/c 2 ) of 
the known D* + -D° (D*°-D°) mass difference. For D* + candidates, the K and it are fit to 
a common vertex, and then the slow tc and D are fit to a second vertex using a beam spot 
constraint. These constraints improve the Am resolution by about 20%. 

To D* candidates we add a lepton candidate. We select electrons using the ratio of 
the energy deposited in the Csl calorimeter to the reconstructed track momentum, the 
shape of the shower in the calorimeter, and the specific ionization in the drift chamber. Our 
candidates lie in the momentum range 0.8 < p e < 2.4 GeV/c. Muon candidates penetrate 
~ 5 interaction lengths. Only muons above about 1.4 GeV/c satisfy this requirement; we 
therefore demand that they lie in the momentum range 1.4 < < 2.4 GeV/c. The charge 
of the lepton must match the charge of the kaon, and, for D* + candidates, be opposite that 
of the slow pion. 

For each candidate we compute 

2E(B)E(D*£) - m| - m(D*£) 2 
C08fe --' = 2|p(B)||p(C-Q| ' (1) 

This quantity helps distinguish signal from D*X£v background and bounds the flight direc- 
tion of the B relative to the D*, which is needed to calculate w. We calculate w for the 
extremes of the B flight direction, and average these two values. The typical resolution in w 
is 0.03. We divide our sample into 10 bins from 1.0 to 1.5, the final bin including candidates 
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up to the kinematic limit of 1.504. For w > 1.25, we suppress background with no loss of 
signal efficiency by restricting the angle between the D* and the lepton. 

At this stage, our sample of candidates contains not only D*lv decays, but also B — > 
D**£i> and non-resonant B — > D*i\tv decays (collectively referred to here as B — > D*Xiv 
decays) and various backgrounds. In order to disentangle the D*iv from the D*Xiv decays, 
we use a binned maximum likelihood fit to the cos &B-D*t distribution. In this fit, the 
normalizations of the various background distributions are fixed, and those for D*iv and 
D*X£u float. 

The distributions of the D*lv and D*Xlv decays are taken from Monte Carlo simula- 
tion p. Radiative B — ► D*£i"y decays, modeled by PHOTOS ||, are treated as signal. 
Non-resonant B — > D*nlv decays are modeled using the results of [^], and B — > D**lv 
decays are modeled using the ISGW2 M form factors. 

We account for five classes of background: continuum, combinatoric, uncorrelated, cor- 
related, and fake lepton. (1) Continuum background from e + e~ — *> qq, which amounts to 
about 3.5% of the candidates in the region —1 < cos0B-D*i < 1 (the "signal region"), is mea- 
sured using off-resonance data taken below the BB threshold. We normalize the cos 9 B -d*i 
distribution to the ratio of on- to off-resonance luminosities, correcting for the small en- 
ergy dependence of the continuum cross section. (2) Combinatoric background events, those 
in which one or more of the particles in the D* candidate does not come from a true D* 
decay, contribute 8% (38%) of the candidates in the D* + (D*°) signal region. We take 
the cos distribution of combinatoric background events from a Am sideband, (155, 

165) MeV/c 2 for D* + and (147, 165) MeV/c 2 for D*°, and normalize using a fit to the Am 
distribution in each w bin. (3) Uncorrelated background, which accounts for approximately 
5% of the candidates in the signal region, arises when the D* and lepton come from the 
decays of different B mesons. Most of this background consists of a D* meson combined 
with a secondary lepton (i.e., from b — > c — > s£u) because primary leptons from the other B 
have the wrong charge correlation in the absence of B°-B° mixing or D* production from 
the hadronization of the c in b — > ccs. We obtain the uncorrelated cos9b-d*i distribution 
from Monte Carlo simulation, normalizing to the inclusive D* production rate observed in 
our data in low and high D* momentum bins, the measured primary and secondary lepton 
decay rates 0, the estimated decay rate for modes in the B — > D*D^K^ family 



and the measured B° — B° mixing rate [IT]. (4) Correlated background events are those in 
which the D* and lepton are daughters of the same B, but the decay was not B — > D*iv or 
B — > D*X£u. The most common sources are B — > D*rv followed by leptonic r decay and 
B — > D*D S followed by semileptonic decay of the D s . This background accounts for fewer 
than 0.5% of the candidates in the signal region and is estimated using Monte Carlo simula- 
tion. (5) Finally, hadrons misidentified as leptons contribute fewer than 0.5% of candidates 
in the signal region. 

Having obtained the distributions in cos 6b-d*z of the signal and background components, 
we fit for the yields of D*iv and D*Xlv decays in each w bin. Representative fits are shown in 
Fig. 0. The quality of the fits is good, as is agreement between the data and fit distributions 
outside the fitting region. The fit results also accurately predict the D* energy distribution 
and the lepton momentum spectrum of the data. 

Given the measured D*£i> yields in ten bins of w, we fit for the partial rate 



5 



60 



m 
o 

c 

> 
HI 



.1 I I I I I I I I I I I I I I I I I I I I I I I I I I I 

(a) 

•,■ Data 

□ D*£v 



:■ D*X^v 

Combinatoric 

□ Continuum 

□ Uncorrelated 

□ Correlated 
2o\- m Fake Leptons 



f\ 




I. ♦ 



Qn 3071201-048 

OU I i i i i | i i i i | i i i i | i i i i | i i i i | i i i i 



60 



40 - 



20 




-10 -8 -6 -4 -2 2 4 -10.0 -7.5 -5.0 -2.5 2.5 5.0 



cos 0, 



B-D £ 



FIG. 1. The candidate yields for 1.10 < w < 1.15 bin with the results of the fit superimposed 
for (a) D* + £~D and for (b) D*°£~i>. The fit uses the region between the arrows. 
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(2) 



where K.(w) is a known function of kinematic variables and T{w) is the form-factor. For our 



fit we use a form-factor parameterization |12| informed by HQET and dispersion relation 
constraints |13|]. We take the form-factor ratios -Ri(l) and -R2(l) from a previous measure- 
ment ]14| that agrees with theoretical expectations |15| . The slope p 2 of the form factor at 
w = 1 is the only shape parameter; it varies in our fit. 

We fit our D*£v yields as a function of w for |V^fe|jF(l) and p 2 . We minimize the sum 



of x for D* + £~u and D 



*0/ 



V. 



each of which compares the measured and expected yields in 



the ten reconstructed w bins. For each mode, 



x = ^ 

i=l 



T 2 

N° 



(3) 



where N° bs is the yield in the i th w bin, Nj is the number of decays in the j th w bin, and the 
matrix accounts for the reconstruction efficiency and the smearing in w. In the above, 



Nj = 4fN r{ 4 S) B D *B D0 T B I dw 



dY 

dw' 



(4) 



where tb is the B 
is the D° —> K 



or B° lifetime |Tl[], Bp* is the D* — > D°tt branching fraction [O, B^o 



7T 



branching fraction 



N 



T(4S) 



is the number of T(4S') events in the 



sample, and / represents /oo or /+-, the T(45) — > B B or B + B~ branching fraction, as 
appropriate. We use the result of |[7| for (/+-//oo)( r +/ r o) as a constraint in the fit, assuming 
foo + /+- = 1- We assume that B~ — > D*°£~u and B° — > D* + £~u have identical partial 
widths and form factors. 
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FIG. 2. The results of the fit to the w distribution. The top figure shows the D* + d~v yields 
(solid circles) with the results of the fit superimposed (histogram). The middle figure shows the 
same for D* £~u. The bottom fi gure displays | V^f,|^ r (?i?) , where the solid circles (squares) are 
derived from the D* + l~v (D* a l~v) yields after correcting for efficiency, smearing, and all terms in 
the differential decay rate apart from | V^fo|j?-"(w) . The curve shows the result of the fit. 



The result of the fit is shown in Fig. |. We find |K fe |.F(l) = 0.0431 ± 0.0013 ± 0.0018, 

p 2 = 1.61 ± 0.09 ± 0.21, and /_) = 0.521 ± 0.012, where the errors are statistical and 

systematic. The fit \ 2 is 16.8/18 d.o.f., and the correlation coefficient between |Vd,|.F(l) and 
p 2 is 0.86, which becomes 0.22 after including systematic correlations. Integrated over w, 
these parameters give T = 0.0394 ± 0.0012 ± 0.0026 ps _1 , implying the branching fractions 
B(B° -> D*+£~P) = (6.09 ± 0.19 ± 0.40)% and B(B~ -> D*°£-p) = (6.50 ± 0.20 ± 0.43)%. 
Separate D* + £rv and D*°£~u fits give consistent results. When fit using the same form 
factor, the present data are consistent with the previous CLEO result [TB|, which analyzed 
a subset of the data reported in this Letter. The results in this Letter supersede our earlier 
measurement. 

The systematic uncertainties are summarized in Table |. The dominant systematic un- 
certainties arise from our background estimations and from our knowledge of the slow-pion 
reconstruction efficiency. 

We test that the combinatoric candidates in the Am sideband have the same cos6b-d*£ 
distribution as the background events in the peak region by applying our procedure to Monte 
Carlo simulated events. The systematic uncertainty is the difference between the results 



7 



obtained using the true Monte-Carlo background and those obtained using the sideband 
subtraction. We also vary the function used to normalize the Am sideband. The main 
source of uncertainty from the uncorrelated background is the branching fraction of the 
B -> D*D^K^ decays, which we vary by ±50%. Smaller effects arise from the primary 
and secondary lepton rates and from the uncertainty in B° — B° mixing. We assess the 
uncertainty arising from the correlated background by varying by 50% the branching fractions 
of the contributing modes. 

A major source of uncertainty for the analysis is the reconstruction efficiency of the slow 
pion from the D* decay. Due to the small energy release in D* decay, the slow-pion momen- 
tum is correlated with w. Thus, for D* + , the efficiency is small near w — 1 for low-momentum 
pions and increases rapidly over the next few w bins, while the efficiency for 7r° — > 77 is more 
uniform in w. We have explored the detection efficiencies as a function of event environment 
(nearby tracks or showers), drift chamber performance (single measurement resolution and 
efficiency), vertexing requirements, calorimeter simulation (noise, nonlinearities, and shower 
simulation threshold), and description of the detector material in our simulation. We vary 
the amount of noise in the calorimeter simulation and introduce possible residual nonlineari- 
ties in the energy scale. These variations are constrained by ^(77) and lateral shower shape 
distributions from an independent sample of tt° candidates from our data. The uncertainty 
in |V^,|J-"(1) is dominated by uncertainties in the number of interaction lengths in the inner 
detector (1.3%) and the vertexing efficiency (1.5%). 

We determine the tracking efficiency uncertainties for the lepton and the K and 7r forming 
the D° in the same study used for the slow pion from the D* + decay. These uncertainties 
are confirmed in a study of 1-prong versus 3-prong r decays. The efficiency for identifying 
electrons (muons) has been evaluated using radiative Bhabha (//-pair) events embedded in 
hadronic events, and has an uncertainty of 2.6% (1.6%). Separate electron and muon analyses 
of our data give consistent results. 

Finally, our analysis requires knowledge of the cos9 B _ D * e distribution of the D*X£u con- 
tribution. This distribution in turn depends on both the poorly known branching fractions 
of contributing modes and their form factors. We note that the B — > D*n£i> and B —>■ D\lv 
modes have the largest and smallest mean cos9b-d*i- We therefore repeat the analysis, first 
using pure B — > D*tt£v to describe our D*X£v decays and then using pure B — > D\£v to 
describe these decays, and we take the larger of the two excursions as our systematic error. 

The form factor ratios R\ and R2 affect the lepton spectrum and therefore the fraction 
of candidates satisfying our lepton momentum requirements. To assess this effect, we vary 
R\ and R2 within their measurement errors, taking into account their correlation. 

Using a recent lattice calculation Jl9| that yields .F(l) = 0.9191q;o35, our result for 
1^1^(1) implies 

\V cb \ = 0.0469 ± 0.0014(stat.) ± 0.0020(syst.) ± 0.0018(theor.). (5) 

Our result is the most precise to date and is somewhat higher than but marginally consistent 
with previous measurements pOfl . However, we note that our ability to reconstruct cos9 B _ D *i 



makes our analysis approximately four times less sensitive to the poorly known D*X£v 
background and allows us to constrain it with the data. This value of | V^^ | is also somewhat 
higher than that obtained from inclusive semileptonic B decays [^TJ. If confirmed, this 



discrepancy could signal a violation of quark-hadron duality. A larger value of \V c b\ affects 
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TABLE I. The fractional systematic uncertainties. 
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constraints on the CKM unitarity triangle, reducing expectations for indirect CP violation 
in the B system. 
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U.S. Department of Energy, the Research Corporation, and the Texas Advanced Research 
Program. 



9 



REFERENCES 



N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. 
Theor. Phys 49, 652, (1973). 

Described more fully in CLEO Collaboration, CLNS 01/1774, to be submitted to Phys. 
Rev. D. 

CLEO Collaboration, Y. Kubota et ai, Nucl. Instrum. Methods Phys. Res. Sect. A 320, 
66 (1992). 

R. Barlow and C. Beeston, Comput. Phys. Commun. 77, 219 (1993). 
R. Brun et al, GEANT 3.15, CERN DD/EE/84-1 (1987). 
E. Barberio and Z. Was, Comput. Phys. Commun. 79, 291 (1994). 
J. L. Goity and W. Roberts, Phys. Rev. D 51, 3459 (1995). 

D. Scora and N. Isgur, Phys. Rev. D 52, 2783 (1995); N. Isgur et al, Phys. Rev. D 39, 
799 (1989). 

CLEO Collaboration, B. Barish et al, Phys. Rev. Lett. 76, 1570 (1996). 
CLEO Collaboration, CONF 97-26 (1997). 

Particle Data Group, D. E. Groom et al, Eur. Phys. J. C 15, 1 (2000). 

I. Caprini, L. Lellouch, and M. Neubert, Nucl. Phys. B 530, 153 (1998) flhep 

ph/9712417|l . 



C. G. Boyd, B. Grinstein, and R. F. Lebed Phys. Rev. D 56, 6895 (1997) |Tep= 
ph/9705252j| . 



CLEO Collaboration, J. Duboscq et al, Phys. Rev. Lett. 76, 3898 (1996). 
M. Neubert, Physics Reports, 245, 259 (1994). 

We use the average of ALEPH Collaboration, R. Barate et al, Phys. Lett. B 403, 367 
(1997) and CLEO Collaboration, M. Artuso et al, Phys. Rev. Lett 80, 3193 (1998). 
These measurements indicate their treatment of final-state radiation. 
CLEO Collaboration, J. P. Alexander et al, Phys. Rev. Lett. 86, 2737 (2001). 
CLEO Collaboration, B. Barish et al, Phys. Rev. D 51, 1014 (1995). 
S. Hashimoto et al, FERMILAB-PUB-01/317-T, |hep-ph/0 110233] . We apply a QED 
correction of +0.007. 

[20] ALEPH Collaboration, D. Buskulic et a/.,Phys. Lett. B 395, 373 (1997); OPAL Col- 
laboration, G. Abbiendi et al, Phys. Lett. B 482, 15 (2000); DELPHI Collaboration, 
P. Abreu et al, Phys. Lett. B 510, 55 (2001); BELLE Collaboration, K. Abe et al, 
Phys. Lett. B 526, 247 (2002). 

[21] CLEO Collaboration, D. Cronin-Hennessy et al, Phys. Rev. Lett. 87, 251808 (2001) 
||hep-ex/0108033 . 



10 



